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Abstract 



. We use a AHQ shell-model calculation with a two-body effective interaction 

OO I derived microscopically from the Reid93 potential to calculate the isospin- 

mixing correction for the ^'^C— ^-^''B superallowed Fermi transition. The effec- 
tive interaction takes into account the Coulomb potential as well as the charge 
dependence of T = 1 partial waves. Our results suggest the isospin-mixing 



■ correction 5c ~ 0.1%, which is compatible with previous calculations. The 



correction obtained in those calculations, performed in a Ohil space, was dom- 



I inated by deviation from unity of the radial overlap between the converted 

' proton and the corresponding neutron. In the present calculation this effect 

is accommodated by the large model space. The obtained 5c correction is 
about a factor of four too small to obtain unitarity of the Cabibbo-Kobayashi- 
■ ^ , Maskawa matrix with the present experimental data. 
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^ . I. INTRODUCTION 

Superallowed Fermi jS transitions in nuclei, {J'" = 0^,T = 1) — > [J'^ = 0"^,T = 1), 
provide an excellent laboratory for precise tests of the properties of the electroweak inter- 
action, and have been the subject of intense study for several decades (cf. Refs. |IHl3 



According to the conserved- vector-current (CVC) hypothesis, for pure Fermi transitions the 
product of the partial half-life, t, and the statistical phase-space factor, /, should be nucleus 
independent and given by 



K 
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where K/{ncf = 27i^ hi2h/{meC^f = 8.120270(12) x 10"^ GeV'S, Gy is the vector couphng 
constant for nuclear (3 decay, and Mp is the Fermi matrix element, Mp = {ipf \ T± \ ipi). By 
comparing the decay rates for muon and nuclear Fermi (3 decay, the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix element |^ between u and d quarks {vud) can be determined 
and a precise test of the unitarity condition of the CKM matrix under the assumption of 
the three-generation standard model is possible 

For tests of the standard model, two nucleus- dependent corrections must be applied 
to experimental ft values. The first is a series of radiative corrections to the statistical 
phase-space factor embodied in the factors 5r and A/j, giving 



/r = /(1 + 5r + A^), (2) 

where 5^ is due to standard, electromagnetic ("inner") radiative corrections (cf. p. 45 in 
Ref. [0) and is what has been referred to as the "outer" radiative correction (cf. p. 47 
of Ref. 0) and includes axial- vector interference terms [§,|l3]. The second correction, which 
is the subject of this work, arises because of the presence of isospin-nonconserving (INC) 
forces (predominantly Coulomb) in nuclei that lead to a renormalization of the Fermi matrix 
element. This correction is denoted by 5c PI,|3|,[T2[| and modifies the Fermi matrix element 
by I Mp P=| MpQ |2 (1 - 5c), where MpQ = [T{T + 1) - TzTz^Y'^ is the value of the matrix 
element under the assumption of pure isospin symmetry. 

With the corrections 5r, A/j, and 5c, a "nucleus-independent" jFt can be defined by 

n = ft{l + 5R + AR){l~5c), (3) 

and the CKM matrix element Vud is given by [|l^ 

2_ 7r^ln2 _ 2984.38(6) s 

where the Fermi coupling constant, Gp, is obtained from muon /3-decay, and includes ra- 
diative corrections. Currently, ft values for nine superallowed transitions have been mea- 
sured with an experimental precision of 0.2% or better P,p^. With these precise mea- 
surements and reliable estimates for the corrections, the CVC hypothesis can be confirmed 
by checking the constancy of the J-'t values for each nucleus, while the unitarity condi- 
tion of the CKM matrix is tested by comparing the average value of Vud with the values 
determined for v^s = 0.2199(17) fig and Vub < 0.0075 (90% confidence level) i.e., 
v"" = vla + vl^ + vlf, = l. 

In the past, the nuclear structure correction 5c has been computed within the framework 
of the nuclear shell model 0,|3,11-13|. In general, the isospin-nonconserving components of 



the nuclear Hamiltonian are small, and can be treated perturbatively. Due to computational 
limitations and uncertainties associated with determining an effective Hamiltonian, almost 
all calculations for nuclei with A> 10 have been performed within a single major oscillator 
shell, e.g., for ^°C the model space spanned by the 0p3/2 and Opi/2 oribitals (;?-shell). Within 
this context, two types of isospin mixing must be accounted for. The first is due to the 
mixing between states that lie within the shell-model configuration space. For example, 
for A = 10, there are 2, 7, and 1 ;?-shell configurations leading to J'^ = 0^ and T = 0, 
1, and 2, respectively. Because of its two-body nature, the INC interaction is composed 
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of isospin operators of rank zero (isoscalar), one (isovector), and two (isotensor), and in 
the case of A = 10 it is capable of mixing together all J'^ = 0^ states. Traditionally, the 
configuration mixing correction is denoted as 6im and in Ref. |TT| it was shown that the 
best estimates for 6im are obtained using an INC interaction that correctly describes the 
Coulomb energy splittings of the binding energies between members of the isospin multiplet, 
e.g., the J'" = 0"*", T = 1 states in ^°C, ^°B, and ^°Be. Of the two types of mixing, 6im is the 
smallest with a magnitude of approximately 0.04-0.1%. 

In addition to the mixing between states contained within the shell-model configuration 
space, mixing with states that lie outside the model space must also be accounted for. In 
particular, the Coulomb interaction can strongly mix one particle-one hole (Ip — Ih) 2hfl 
excitations, e.g., Ops/i 1^3/2, into the ground state. In previous works, excitations of this 
type were accounted for by examining differences in the single-particle radial wave functions. 
Indeed, for closed-shell configurations, mixing with Ip — Ih states is properly accounted for 
at the level of Hartree-Fock. Hence, the second correction to the Fermi matrix element, 
denoted by Sjio, was estimated by evaluating the mismatch in the radial overlap between 
the single-particle wave functions of the converted proton and the corresponding neutron. 
The explicit details for the calculation of 6ro, which involve a sum over intermediate A — 1 
parent states that then determine the proton and neutron separation energy for the radial 



wave function, are given in Refs. PJTl[|. For the most part, 6^0 is found to be the larger of 
the two components (with 6c = 5ro + hu) and has a magnitude of the order 0.1-0.8%. 

At present, two methods for evaluating 5ro are espoused. The first (THH) [0 uses 
Woods-Saxon (WS) radial wave functions, while in the second (OB) [pI PIpTSf , Hartree-Fock 
(HF) wave functions are employed. Generally speaking, the two methods yield approxi- 
mately the same dependence on nucleon number A, but the HF values are systematically 
smaller by 0.1% for the magnitude of the correction. The reason for the difference lies in the 
HF mean field. The principal effect of the Coulomb interaction is to push the proton wave 
functions out relative to the neutrons, hence, providing a mismatch in the radial overlap. In 
Hartree-Fock, however, the proton and neutron mean fields are coupled, and the Coulomb 
interaction actually induces an attractive isovector mean field between the protons and neu- 
trons. In effect, the Coulomb interaction pushes the protons out, but because of the strong 
interaction, the protons pull the neutrons out with them, hence, reducing the magnitude of 
the radial overlap mismatch. 

When all known corrections, i.e., 5r, A/j, 5roi and (5/m, are applied to the nine experi- 



mental data 0,0, it is found that the values are essentially constant within the limits 
of uncertainty but the unitarity limit is violated at the level of approximately 0.4(1)% or 
0.3(1)% for the OB and THH corrections, respectively. In addition, preliminary data from 
a new experiment for ^°C [0 leads to an jFt value that is significantly smaller than that 
of Ref. fl^, and has been interpreted as possible evidence for an, as yet, unaccounted for 



correction that might lead to satisfying the unitarity condition of the CKM matrix. In ad- 
dition, it must be admitted that the present separation between the configuration mixing 
and radial overlap contributions to 5c is somewhat unsatisfying. A much better approach 
would be to perform a shell-model calculation that includes several hVL excitations, so that 
both corrections would evaluated on the same footing and simultaneously. Because of re- 
cent improvements in computational capabilities and the ability to determine an effective 
model-space Hamiltonian based on realistic nucleon-nucleon interactions, it is now possible 
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to perform such a calculation for the lightest of the nine accurately measured transitions. 
We report here the results of large-basis shell-model calculations that include excitations up 
to 4hQ for A = 10 nuclides, with an emphasis on evaluating the isospin-mixing corrections 
to the matrix element for the Fermi decay of ^"C 

The organization for the paper is as follows. First, in Section |1| we discuss the shell- model 
Hamiltonian with a bound center-of-mass, the method used to derive the starting-energy- 
independent effective interaction, and the renormalization of the transfer operator. Results 
of the Fermi matrix element calculations are presented in section |ITl| and concluding remarks 
are given in section |V]. 



II. SHELL-MODEL HAMILTONIAN AND THE EFFECTIVE INTERACTION 



In our calculation we use the one- plus two-body Hamiltonian for the A-nucleon system. 



I.e., 



i=i i<j 



(5) 



where m is the nucleon mass and VN(^i — rj) the nucleon-nucleon interaction, modified by 
adding the center-of-mass harmonic-oscillator potential ^Amil'^R^, R = \ j2t=i'^i- This 
potential does not influence intrinsic properties of the many-body system. It provides, 
however, a mean field felt by each nucleon and allows us to work with a convenient harmonic- 
oscillator basis. The modified Hamiltonian, depending on the harmonic-oscillator frequency 
Q, may be cast into the form 



E 



if 1 

2m 2 



A 

E 

Kj L 



(6) 



The one-body term of the Hamiltonian ([ 



X]j=l Pi) 



is then re-written as a sum of the center-of-mass 
and a term depending on relative coordinates 



term H^^ — , 

only. Shell-model calculations are carried out in a model space defined by a projector P. 
In the present work, we will always use a complete NhVL model space. The complementary 
space to the model space is defined by the projector Q = 1 — P . In addition, from among 
the eigenstates of the Hamiltonian (|^), it is necessary to choose only those corresponding 
to the same center-of-mass energy. This can be achieved by projecting the center-of-mass 
eigenstates with energies greater than upwards in the energy spectrum. The shell-model 
Hamiltonian, used in the actual calculations, takes the form 



H 



P0 



A 

E p 

i<j=i 



2Am 



+ 



2T 



r,- - r,- 



P + Y.P 

i<j 



~ ~2A 



P 



J eff 



(7) 



where /3 is a sufficiently large positive parameter. 

The effective interaction introduced in Eq. (|^) should, in principle, exactly reproduce 
the full-space results in the model space for some subset of states. In practice, the effective 
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interactions can never be calculated exactly as, in general, for an A-nucleon system an A- 
body effective interaction is required. Consequently, large model spaces are desirable when 
only an approximate effective interaction is used. In that case, the calculation should be less 
affected by any imprecision of the effective interaction. The same is true for the evaluation 
of any observable characterized by an operator. In the model space, renormalized effective 
operators are required. The larger the model space, the less renormalization is needed. 

Usually, the effective Hamiltonian is approximated by a two-body effective interac- 
tion determined from a two-nucleon system. In this study, we use the procedure as de- 
scribed in Ref. |]T7[. To construct the effective interaction we employ the Lee-Suzuki 
p!8| similarity transformation method, which gives an interaction in the form ^2^5-^2 = 
P2VP2 + P2VQ2U!P2, with u the transformation operator satisfying u = Q20jP2- The projec- 
tion operators P21Q2 = 1 — P2 project on the two-nucleon model and complementary space, 
respectively. Our calculations start with exact solutions of the Hamiltonian 



^ Hg + V,"" = + + rf) + V{n - ^2) - - rl)^ . (8) 



which is the shell-model Hamiltonian (^) applied to a two-nucleon system. We construct the 
effective interaction directly from these solutions. Let us denote the relative-coordinate two- 
nucleon harmonic-oscillator states, which form the model space, as |ap), and those which 
belong to the Q-space, as \q.q). Then the Q-space components of the eigenvector \k) of the 
Hamiltonian (H) can be expressed as a combination of the P-space components with the help 
of the operator oj 

{aQ\k) = J2{(^QH(^p){(^p\k) ■ (9) 

ap 

If the dimension of the model space is dp, we may choose a set /C of dp eigenevectors, 
for which the relation (|^) will be satisfied. Under the condition that the dp x dp matrix 
{ap\k) for \k) G /C is invertible, the operator u can be determined from (|^). In the present 
application we select the lowest states obtained in each channel. Once the operator u is 
determined the effective Hamiltonian can be constructed as follows 



{'lp\H2es\ap) = 



{-fp\k)Ek{k\ap) + J2{^p\k)Ek{k\aQ){aQ\uj\ap) 



(10) 



This Hamiltonian, when diagonalized in a model-space basis, reproduces exactly the set /C 
of dp eigenvalues E/^. Note that the effective Hamiltonian is, in general, quasi-Hermitian. 
It can be hermitized by a similarity transformation determined from the metric operator 
P2(l + uj'^uj)P2. The Hermitian Hamiltonian is then given by 



2eff 



1/2 



P2{1+UJ^UJ)P2 H2,S P2{1+UJ^UJ)P2 . (11) 



-1/2 



Finally, the two-body effective interaction used in the present calculations is determined 
from the two-nucleon effective Hamiltonian (|ll]) as V^s = -f^2cfT ~ -f^02- Note that we dis- 



tinguish the two-nucleon system projection operators P2,Q2 from the A-nucleon system 
operators P, Q. 
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To at least partially take into account the many-body effects neglected when using only 
a two-body effective interaction, we employ the recently introduced multi-valued effective 
interaction approach ||2^. As a consequence, different effective interactions are used for 



different excitations. The effective interactions then carry an additional index indicating 
the sum of the oscillator quanta for the spectators, A^sps, defined by 

N = N - N - N ■ = N' - N - N ■ (12) 

-"'sps -''sum -''a -''spsmm -''sum 7 -''spsmm ; \ ) 

where Agum and Ag'um total oscillator quanta in the initial and final many-body 

states, respectively, and A^, and A^ are the total oscillator quanta in the initial and final 
two-nucleon states \a.) and I7), respectively. A'spsmin is the minimal value of the spectator 
harmonic-oscillator quanta for a given system. Here, for A=10, Aspsmin = 4. Different sets 
of the effective interaction are determined for different model spaces characterized by A'spg 
and defined by projection operators 

P2(A^sps) = l-Q2(A^sps) . (13b) 

In Eqs. (0), A^max characterizes the two-nucleon model space. It is an input parameter 
chosen in relation to the size of the many-nucleon model space. This multi-valued effective- 
interaction approach is superior to the traditional single-valued effective interaction, as con- 



firmed also in a model calculation |2l|] 



Our goal in this study is to evaluate the Fermi matrix element 

Mp = (^°B, 0+l|T_|i°C, 0+1) , (14) 

which is equal to \/2 for an isospin-invariant system. Note that for a system with isospin 
breaking, the isospin-lowering operator T_ should be renormalized in a similar way, as 
the interaction used for calculation of the eigenstates appearing in Eq. (|l^) . In fact, we 



can apply the formalism described in Refs. |22| to construct a two-body effective operator 



(T_)efr consistent with the two-body effective interaction derived above and exact for the 
two-nucleon system. Then we could use such an operator in the A-body calculation. We 



studied such a possibility in a solvable-model calculation as described in Refs. Here, we 
did two-nucleon calculations with the effective T_ operator. The observed renormalization 
of the bare operator for the model spaces of the size used in our calculations, was, however, 
insignificant compared to the other effects as described further. Therefore, in the A-body 
calculations we used the bare T_ operator. 



III. APPLICATION TO THE A=10 SYSTEM WITH ISOSPIN BREAKING 

In order to evaluate the Fermi matrix element [Eq. (0)], we apply the formalism outlined 
in section || for A=10 nuclei. In the calculations we use the Reid93 nucleon-nucleon potential 
p3| and consider the following isospin-breaking contributions. First, the Reid93 potential 
differs in the T = 1 channels for pn and pp, nn systems, respectively. Second, we add the 
Coulomb potential to the pp Reid93 potential. Consequently, using the Eqs. (|9D- (|ril) , we 
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derive different two-body effective interactions for the pn, pp, and nn systems. No other 
mechanisms for isospin breaking are considered. 

As we derive the effective interaction microscopically from the nucleon-nucleon interac- 
tion, the number of freely adjustable parameters in the calculation is limited. 

First, we have the choice of the model-space size in the shell-model diagonalization. That 
is, however, constrained by computer capabilities. The largest model space we were able 
to use was the space allowing all 4hQ excitations relative to the unperturbed ground state. 
Most of the calculations were done in the m-scheme using the Many-Fermion-Dynamics 
Code extended to allow the use of different pn, pp, nn interactions. We also performed 
some calculations with the OXBASH shell-model code In the m-scheme, the dimensions 
associated with and are 581,740 and 430,137, respectively. To study the dependence 
on the model-space size, we performed calculations in the 2hQ space as well. In that space, 
the dimensions drop to 14,502 and 10,111, respectively. 

Second, we have the choice of the two-nucleon model space used for the evaluation of the 
effective interaction. This is related to the many-nucleon model-space size, and, in principle, 
is determined by that size. Traditionally, however, the Q = space used to determine the 



G-matrix does not, necessarily, coincide with the many-particle model space [p6| , p7| . In 
our calculation, the two-nucleon model space is characterized by a restriction on the number 
harmonic-oscillator quanta Ni < N^a.^, < A^max, (A^i + A^2) < A^max- Here, Aj = 2nj + Zj is 
the harmonic-oscillator quantum number for the nucleon i,i = 1, 2. This type of restriction 
guarantees an orthogonal transformation between the two-particle states and the relative- 
and center-of-mass-coordinate states. For the present AhQ calculation, the choice of A^max = 
6 appears to be appropriate. However, it has been observed in the past p!7|j28|j29[] that when 



the Lee-Suzuki procedure combined with the G-matrix calculation according to Ref. |2^ 



(which is equivalent to the procedure we are using) is applied to calculate the two-body 
effective interaction, the resulting interaction may be too strong. This is, in particular, true, 
when the multi-valued approach is used. Several possible adjustments were discussed to 
deal with this problem [|l3,0 and amounted to introducing an extra parameter. In the 
present calculations, we do not introduce any new parameter, but rather we treat A^max ^ 
a free parameter and use Amax = 8 for the AhQ calculations and A^max = 6 for the 2hQ 
calculations, respectively. With this choice we obtain quite reasonable binding energies for 
the studied nuclei. We have also performed several 2hQ test calculations with single-valued 
interactions that were derived following Ref. 0], as opposed to the multi- valued interaction 
discussed in the previous section. To obtain reasonable binding energies with the single- 
valued interaction we do not have to change the A^ax value from that corresponding to the 
many-nucleon space, e.g., A^max = 4 for the 2hQ calculation. This difference of treatment 
of the two types of interactions follows from the fact that the overall strength of the single- 
valued interaction is weaker. 

Third, our results depend on the harmonic-oscillator frequency Q. We have studied this 
dependence by performing calculations for the values HQ = 14, 15.5, and 17 MeV. 

Let us also mention one important feature of the present approach. For both the multi- 
valued and the single- valued interactions our calculations do not break the separation of the 
center-of-mass and the internal relative motion. In particular, a variation of the parameter 
f3 introduced in Eq. (^) does not change the eigenenergies and other characteristic of the 
physical states. This is so due to the choice of a complete Nhfl many-nucleon space and 
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the triangular two-nucleon model space for deriving the effective interaction as well as the 
due to the procedure used to derive the effective interaction. 

In Figs. |I], 0, and ^ we present the experimental and calculated spectra of for 
hQ = 14, 15.5, and 17 MeV, respectively, for the 2hQ and AhQ model spaces. In general, we 
observe an overall improvement in the spectra with the enlargement of the model space in 
all three cases. Also, the 4hQ calculations exhibit more stability with regard to changes in 
the harmonic-oscillator frequency than do the 2hQ results. The agreement with experiment 
improves when going from hQ = 14 MeV to HQ = 17 MeV; in particular for the ground 
state and the lowest states. In fact, from Fig. ^ we find that a very reasonable description 
of the spectra is obtained for hQ = 17 MeV. 

In Table |[ the overall behavior with respect to hQ is illustrated. In general, we observe 
a reasonable reproduction of the binding energy, with a moderate decrease occuring for 
increasing hQ. Using free-nucleon effective charges, we find that although the quadrupole 
moment for the 3^0 state is underestimated considerably, the magnetic dipole moment is 
well reproduced. In addtion, the point-proton rms radius exhibits a fairly strong dependence, 
and increases with decreasing HQ. For the rms radius, we find that the best agreement with 



experiment |^2[ is achieved for hQ = 14 MeV. 

From the point of view of the beta decay of ^°C, a good description of the T = 1 states 
is important. From Figs. |l|-^ we can see that the calculated ^'^B T = 1 states have the 
right relative positions and are reasonably stable with variations of both the model-space 
size and hQ. We have also performed AhQ calculations for ^°Be to study the splitting of 
the isospin analog states in the whole isospin-multiplet - ^°B - ^°Be. The experimental 
ground state splitting between and ^°Be is 4.66 MeV, while our calculated values are 
4.68, 4.83, and 4.94 MeV for hQ = 14, 15.5, and 17 MeV, respectively. The best agreement 
with experiment is achieved for hQ = 14 MeV, where the calculated rms point-proton radius 
is also in agreement with the experimental value. On the other hand, the splitting between 
the 0"'"1 states of and ^°B, which is experimentally 2.69 MeV, is overestimated in our 
calculations by 8, 11, and 14% for the hfl = 14, 15.5, and 17 MeV calculations, respectively. 
Since the correct - ^°Be splitting is obtained for hQ = 14 MeV, the excess in the ^"C - 
^°B splitting suggests that the isospin breaking due the strong T = 1 force may be too large. 
One possible explanation is that our approach for deriving the effective interaction tends to 
exaggerate the differences between the pn and nn, pp potentials. Such an artificial effect 
should decrease with increasing model-space size. On the other hand, it is also possible that 
the Reid93 potential itself overestimates differences between the pn and pp, nn systems in 
the T = 1 channel. For the most part, we find the best overall agreement for the rms point 
proton radius, binding energy, and Coulomb energy splitting for hQ = 14 MeV. Given that 
the isospin mixing is largely driven by the Coulomb interaction, which is then dependent on 
the size of the nucleus, we feel that the best value for the isospin- mixing correction to the 
Fermi matrix element will be achieved for HQ = 14 MeV. 

The most important results of our study are also summarized in Table | in the last two 
lines. The calculated isospin- mixing corrections Sq = I — ^^J^ , in %, are presented for all 
three choices of Q and for both 4hQ and 2hfl model spaces. Again, a correlation between the 
radius and the isospin-mixing correction is clearly observed, as Sc decreases with increasing 
radius. This is simply understood in terms of a larger radius implying weaker Coulomb 
effects. On the other hand, with an increase in the model-space size, a significant increase 
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in the isospin-mixing correction is apparent. This is due to the fact that in the larger model 
space, the excitation energies of the Ip-lh O"*"! states decrease, hence leading to greater 
mixing. For this reason, the more realistic multi- valued effective interaction is important. 
We have also performed test calculations with the single- valued interaction in the 2hfl space, 
and found 6c to be smaller by approximately 30%. 

Our 4hQ results suggest an isospin-mixing correction 5c ~ 0.08—0.1%. This is compatible 
with the previously published value of 6c ~ 0.15(9)% by Ormand and Brown []T^. That 
value was a sum of two contributions. First, about 0.04% came from the shell-model wave- 
function renormalization due to the isospin mixing and was obtained in a Ohfl shell-model 
calculation using phenomenological effective interactions. Second, the amount 0.09% was 
due to the deviation from unity of the radial overlap between the converted proton and the 
corresponding neutron. This effect was attributed to the influence of states lying outside 
the OhQ space. The radial wave-functions were obtained in a Hartree-Fock calculation using 
Skyrme-type interactions. Because we use a multi-configuration model space in the present 
calculation, we should have both effects included consistently at the same time. 

Another important factor in the calculation is the position of the 2hQ states. As dis- 
cussed before, the position of the Ip-lh states influences the ground-state isospin mixing. 
Unfortunately, the excitation energy of these states is not known experimentally. However, 
in our calculations the multi-valued effective interaction is used and a more realistic de- 
scription of these states should be obtained; especially, in the ihQ model space. On the 
other hand, in an analogous calculation for ^He, it was observed that an 8hQ model space 
is needed to get the 2tiQ dominated 0^ state close to the experimental excitation energy 
I3,|20|. There are states like 1+0 at 5.18 MeV in ^°B or 0+1 at 6.18 MeV in ^°Be, which are 



believed to be 2p-2h, 2hil excitations. We do not observe any such states below 7.5 and 12 
MeV, respectively, in our calculations. The first excited 0+1 state in ^''Be obtained in the 
AhQ calculation with hQ = 14 MeV lies at 9.8 MeV. It is, however, predominately a Ohfl 
state. There can be two reasons why we do not get such states. First, these states have not 
yet converged in the Lanczos procedure. Second, and more likely, the 4hQ model space is too 
small for the right description of the 2hQ excitation states. Therefore, it would be desirable 
to extend the present calculations to a larger, e.g., 6hQ model space. Unfortunately, due to 
the computational limitations, it is not possible at this time to perform a calculation of this 
magnitude. However, from Table | we observe an increase of 6c by ~ 0.03% between the 
2hQ and the 4hQ calculation. Therefore, we might expect an increase of similar magnitude 
for an increase of the model-space size beyond 4hQ. Therefore, the more realistic value 
of the isospin-mixing correction from our calculation would be 6c ~ 0.12(3)%, where the 
uncertainty is estimated from the change in 6c obtained when using an increased model 
space. 



IV. CONCLUSIONS 

The effects of isospin mixing on the transition matrix element for the superallowed Fermi 
/9-decay of were estimated within the context of a large-basis, shell-model calculation. 
The calculations were performed assuming no closed core and an effective interaction based 
on a realistic two-body nucleon-nucleon interaction, while including the Coulomb interaction 
between protons. Contrary to previous estimates for the isospin corrections, this calculation 
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was carried out within a model space that inchidcd many hfl excitations. As a consequence, 
the conventional configuration mixing and radial mismatch contributions were evaluated 
within a unified framework, simultaneously and the usual separation was not necessary. With 
regard to parameters used within the calculation, we find a correlation between the isospin- 
mixing correction and the Coulomb splitting between the isotopic multiplets, which, in turn, 
is governed by the nuclear size through the oscillator parameter. Given that the isospin- 
mixing correction is primarily a Coulomb effect, the best value for Sc is taken to coincide 
with the oscillator parameter that correctly reproduces the Coulomb splittings. With regard 
to the model-space size, a clear improvement (or an indication towards convergence) in most 
observables is evident when the size of the model space is increased from 2hfl to AhQ, but 
Sc is found to increase by only 0.03% (in magnitude) in this case. Hence, our final estimate 
for 5c is taken to be 0.12(3)% (where the 4hQ result has been increased by 0.03% to account 
for the possible effects of an increased model space). This result also happens to be in 
excellent agreement with the previous estimates that relied on the conventional separation 
of the configuration and radial mismatch contributions. Finally, we note that the magnitude 
of the isospin-mixing correction obtained in our calculation does not lead to a resolution to 
the deviation from unitarity for the Cabibbo-Kobayashi-Maskawa matrix. 
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FIGURES 



FIG. 1. The experimental and calculated excitation spectra of B. The results corresponding 
to the model-space sizes of ihO, and 2hO, relative to the ground-state configurations are presented, 
respectively. The harmonic-oscillator energy of hfl = 14 MeV was used. 



FIG. 2. The same as in Fig. || for the harmonic-oscillator energy of ^0 = 15.5 MeV. 
FIG. 3. The same as in Fig. |^ for the harmonic-oscillator energy of hCl = 17 MeV. 
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TABLES 





Exp 


hn = 14 MeV 


nn = 15.5 MeV 


nn = 17 MeV 


^b(IOB) 


64.75 


63.61 


62.78 


61.53 


<9(3+0) 


8.47(6) 


5.85 


5.64 


5.52 


/i(3+0) 


1.80 


1.86 


1.85 


1.85 


/.(l+O) 


0.63(12) 


0.84 


0.84 


0.84 




60.32 


58.68 


58.19 


56.83 




2.31 ±0.03 


2.28 


2.21 


2.17 


Sc{4nn)[%] 




0.084 


0.091 


0.097 


5ci2nn)[%] 




0.055 


0.061 


0.067 



TABLE L Experimental and calculated binding energies, in MeV, magnetic moments, in //n, 
and quadrupole moments, in e fm^, of ^'^B. Also the experimental and calculated binding energies, 
in MeV, and the point proton radius, in fm, of ^''C are presented. The results correspond to the 
4:hfl calculations. In addition, the isospin-mixing correction 5c, in %, is shown as obtained both 
in the 47i$7 calculations and the 27iO calculations. Results of three different calculations with the 
harmonic-oscillator parameter taken to be HQ = 14, 15.5, 17 MeV, respectively, are presented. The 
effective interaction used was derived from the Reid 93 nucleon-nucleon potential. The experimental 



values are taken from Refs. |31,32|. 



14 



T 1 

r 

T 1 

r 

0^ 1 

r 

3^0 



B 



liQ=14McV 



Exp 4tiQ 2\Q. 



T 1 

r 

T 1 

r 

0^ 1 

r 

3^0 



B 



liQ=15.5McV 



Exp 4tiQ 2\Q. 



T 1 

r 

T 1 

r 

0^ 1 

r 

3^0 



B 



hQ=17MeV 



Exp 4tiQ 2\Q. 



